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Phenomena around us seem to evolve only in the 
forward direction of time, despite the fact that the 
underlying microscopic laws of physics are invariant 
under time reversal apart from a certain rare inter- 
action in particle physics This enigma concern- 
ing the arrow of time is naturally associated with the 
second law of thermodynamics, which dictates that 
the entropy of an isolated system never decrease 
The question of whether and, if so, how the second 
law emerges from quantum mechanics has long been 
tackled |3]. However, the current understanding re- 
mains unsatisfactory because most studies deal with 
open systems coupled to the surrounding environ- 
ments where the time evolution of the system alone 
is no longer time-reversal symmetric (j-lfll- Here we 
show that the second law emerges from an isolated 
pure quantum state through its unitary evolution. 
The physics behind the second law is the quantum- 
mechanical energy uncertainty associated with every 
finite-time operation that makes a large number of 
many-body energy eigenstates indistinguishable, re- 
sulting in a non-decreasing entropy of the system 
Our result establishes the quantum-mechanical def- 
inition of the thermodynamic entropy, opening up 
the way to study thermodynamics in isolated quan- 
tum systems, where persistent quantum coherence is 
expected to cause hitherto unexplored effects ^ 01. 
In fact, independently of details of the system and 
processes, quantum coherence between many-body 
eigenstates is shown to bring about a universal many- 
body correction to the entropy. Such new effects are 
expected to be experimentally observed using ultra- 
cold atoms or ions 

The second law of thermodynamics asserts that the 
entropy of an isolated system never decreases during 
any thermodynamic process. On the other hand, quan- 
tum mechanics tells us that an isolated system is de- 
scribed by a pure quantum state and that its time evo- 
lution is unitary and hence reversible. Then, the crucial 
question is: can such a non-decreasing quantity emerge 
from pure quantum states? 

We answer this question in the affirmative by propos- 
ing that the proper quantum-mechanical definition of 
the thermodynamic entropy is given by 
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FIG. 1: Summary of the main result. With the control 
parameter in the Hamiltonian held fixed, the pure quantum 
state reaches equilibrium (blue). Then a thermodynamic pro- 
cess (purple), which consists of the operation (yellow) and 
equilibration (red), causes a transition to a different pure state 
at equilibrium. During the thermodynamic process, the en- 
tropy never decreases, S < S' . Since both the initial and final 
states are pure, repeated application of this argument leads to 
the second law of thermodynamics, S < S' < S" < ■ ■ . 



matrix p in the eigenbasis of the Hamiltonian and Boltz- 
mann's constant is set to unity 111. El]. It quantifies the 
energy uncertainty and correctly describes the thermo- 
dynamic entropy. It is because it reduces to the von 



where pnn^ are the diagonal elements of the density 



Neumann entropy for the canonical ensemble 11311 . We 
note that the von Neumann entropy is inconsistent with 
the second law of thermodynamics because it is invari- 
ant under any unitary evolution in an isolated quantum 
system and vanishes for pure quantum states. 

Before stating our main result, we recall the zeroth 
law of thermodynamics for pure quantum states, which 
states that a pure state reaches equilibrium under uni- 
tary evolution ITsll . This fact has only recently been ver- 
ified theoretically |[l6l,[l7ll and experimentally lisll , and 
the closely related result concerning quantum ergodic- 
ity has been unearthed fl9^. For equilibration to occur, 
the pure state must involve a large number of energy 
eigenstates, which are coherently superposed. 

Our main result is that the entropy ^ never de- 
creases during any thermodynamic process in large 
systems as follows. Let a thermodynamic process be 
composed of the operation and equilibration, causing 
a transition between two pure states at equilibrium as 
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FIG. 2: Geometrical interpretation of the entropy increase, (a) Each density matrix is represented by a point on or inside the 
Bloch sphere, (b) The entropy {T} depends only on the z-component and attains the maximum value of In 2 on the equatorial 
plane, (c) A pure state traces a circular quantum orbit perpendicular to the z-axis under unitary evolution for a fixed Hamil- 
tonian. The projection of this orbit onto the z-axis gives Pdec/ which retains no coherence between |0) and (d) A quantum 
operation and the associated change in the energy eigenbasis, which tilt the axis of the orbit. The entropy of the decohered state 
i^dec iricreases during the process The entropy increases for more than half of the pure states on the orbit. 



illustrated in Fig. [T] A control parameter in the Hamil- 
tonian, e.g. a magnetic field, is varied during the oper- 
ation and kept constant during equilibration. Let S and 
denote the entropies before and after the thermody- 
namic process, respectively. Then, assuming a uniform 
distribution about the timing at w^hich the operation is 
performed, we obtain 

S < (2) 

for large systems with almost unit probability. The 
equality holds for quasistatic, or infinitely-slow, pro- 
cesses. Note that the initial state need not be assumed 
to belong to any special statistical ensemble such as the 
canonical ensemble [|2Q|-|22|i and that the quantum states 
before and after the operation are both pure states. 
Thus, by applying the result (|2]) to a sequence of op- 
erations, we obtain the second law of thermodynamics 
(see Fig.[lJ. 

The physics behind the entropy increase (0 is the 
time-energy uncertainty principle in quantum mechan- 
ics, which requires that any operation during a time At 
causes the energy uncertainty about h/At, where h is 
Planck's constant |7]. In large systems, there exists an 
exponentially large number of many-body eigenstates 
within the energy window h/At which cannot be distin- 
guished and hence contribute to the entropy increase. 

To illustrate the physical origin of the entropy in- 
crease, let us consider a qubit with two energy eigen- 
states |0) and Each quantum state, which is de- 
scribed by a 2 X 2 density matrix p, is represented by 
a point on or inside a unit sphere called the Bloch 
sphere [23] (see Fig. [2^). Here the x-, y-, and z- 
coordinates are tr[p(7fl] {a = 1,2,3), where CaS are the 
Pauli matrices. The entropy of each state depends only 
on its z-component and monotonically increases from 
the poles to the equatorial plane (see Fig.lJJ)). 

During the unitary evolution with the time- 
independent Hamiltonian, a pure state traces a circle 



around the z-axis which we call the quantum orbit. The 
corresponding mixed state pdec is defined by eliminat- 
ing the coherence between |0) and |1) or, equivalently, 
by projecting the quantum orbit onto the z-axis (see 
Fig.Efc). 

A quantum operation, together with the change of 
the energy eigenbasis, causes a unitary transformation 
on the quantum orbit, which is represented by a rota- 
tion on the Bloch sphere. The quantum orbit is trans- 
formed to a circle around a different axis with p^^^ 
changing into p^^^ (see Fig.|2li). 

From this consideration, it is obvious that the entropy 
of jO^ec never decreases after the operation because jO^ec 
approaches the equatorial plane. For the pure states 
on the orbit, the entropy is more likely to increase than 
decrease (see Fig. |2li). The inequality Q thus holds 
with more than 50% probability. We note that this argu- 
ment is valid only if the quantum orbit never crosses the 
equatorial plane. However, this condition, in practice, 
does not impose any restriction in the thermodynamic 
limit because the equatorial plane, which corresponds 
to the infinite temperature, cannot be reached. 

The asymmetry in the entropy change originates 
from the fact that the quantum orbit is restricted to a 
circle around the z-axis, which breaks the invariance 
of the Hilbert space under any unitary transformation. 
No matter how complex the system might be, the en- 
ergy is conserved and the Hamiltonian is well-defined 
and hence gives the preferred basis, which restricts the 
possible quantum orbits. It is remarkable that the en- 
tropy increase thus results from the conservation of the 
energy. 

The absence of Maxwell's demon ||23-[27|l is also es- 
sential for the entropy increase. It is because there ex- 
ist some timings at which the operation decreases the 
entropy even though the probability is less than 50% 
for a qubit system and vanishingly small for a high- 
dimensional Hilbert space as shown below. Thus, un- 
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less one is able to know all the information on the pure 
state and do infinitely-fine tuning on the timing of the 
operation, the entropy increase is unavoidable. 

For the general case, where the dimension of the 
Hilbert space is d, the entropy increase becomes much 
more probable than for a qubit system. The deco- 
hered state jOdec is defined by eliminating all the coher- 
ence between many-body eigenstates involved in the 
pure state. Then the entropy of p^^^ increases during 
any thermodynamic process in which p^^^ changes into 
P'dec' ^ ^ ^dec^ ^here S^j^^ = S{p'^J, because pdec ap- 
proaches the equatorial plane (see section I of Supple- 
mentary Information for details). For the pure states, 
ignoring a non-extensive correction (see Methods Sum- 
mary), we obtain = S^^^ and thus S < S' except 
for the probability of the order of d~^^^. In many-body 
systems, this is exponentially small with the number of 
particles. 

Finally we demonstrate our result for a system of 
hard-core bosons on a lattice. We place five bosons on 
twenty sites, which are arranged in a 4 x 5 rectangle 
(see Fig. (3^), where the boundary is free. The kinetic 
energy is given by H^in = -/E(z,y)(&- + bpi), where 

bj and bj represent the creation and annihilation oper- 
ators of the boson on site i and the sum is taken over 
all the nearest neighbors, which means that each boson 
can hop to adjacent sites and the energy accompanied 
by the hopping is — /(< 0). The hard-core condition is 
imposed by requiring that more than one particle can- 
not reside on any single site. We introduce a linear po- 
tential h along the y-axis, which gives a tilting potential 
energy, Hpot = h J^i yfii, where ni is the number of par- 
ticles on site i and y/ = 0, 1, 2 and 3 is the y-coordinate 
of the site i as illustrated in Fig.jS^. Applying this linear 
potential is equivalent to tilting the rectangle up at one 
edge in the gravitational field as illustrated in Fig.jSJ^. 

Initially, the linear potential is switched off and the 
quantum state is prepared in the 1000th energy eigen- 
state and thus stationary. The entropy of this state is 
S = 0, which is the minimum. The first operation is 
the quantum quench at time f = 0, which suddenly 
switches on the linear potential to h = / as illustrated 
in Fig.|3j:. The entropy after the quench is = 7.15 and 
obviously S < holds. 

The second operation is the quench at f = r, which is 
a sudden switch-off of the linear potential from h = J to 
(see Fig.lSt). The entropy after the second quench S^^ 
is shown in Fig. |3li. For the quench after equilibration, 
S" coincides with S'^^^ apart from a finite-size correc- 
tion and thus S' < S^^ holds, where S'^^^ is the entropy 
after the second quench assuming that all the coherence 
between many-body eigenstates is lost after the first 
quench. It is remarkable that the finite-size correction 
is universal and given by 7 — 1, where 7 = 0.5772 . . . 
is Euler's constant (see Methods Summary), since it is 
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FIG. 3: Numerical demonstration of the main result, (a) 

The model. Five hard-core bosons are placed on the 4x5 
sites whose boundary is free. The bosons can hop to adja- 
cent sites, (b) Applying a linear potential h along the y axis 
tilts the potential energy of the rectangle up at one edge in 
the gravitational field g. (c) The protocol of the linear po- 
tential. For t < 0, we prepare the 1000th energy eigenstate, 
which is stationary. Two operations performed at ^ = and 
T are the switch-on and -off of h. (d) Entropy change after 
the second quench, or the switch-off of h at t = r, plotted 
against the quench time r. The shaded area corresponds to 
the quench performed before equilibration (the correspond- 
ing region is indicated by the red arrow in (c)).. The entropy 
in the unshaded area is independent of r and equal to the 
one obtained for the decohered state except for the universal 
correction 7 — 1. The inset shows a schematic figure of the 
entropy increase. A quantum state is confined in a constant- 
entropy surface after equilibration. 
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independent of the details of the system and processes. 

The inset of Fig. [Sji schematically illustrates a situ- 
ation in which the quantum state is always pure and 
never converges to the decohered state during the uni- 
tary evolution. Nevertheless, the trajectory eventually 
approaches a constant-entropy surface whose entropy 
is the same as that of the decohered state. Our find- 
ing is that the quantum coherence between many-body 
eigenstates gives no more than a universal correction, 
which is negligible in the thermodynamic limit. 

We remark the difference between equilibration and 
decoherence between the two quenches. So long as 
quantum-mechanical observables are concerned, we 
cannot distinguish them after the second quench flil, 
[itI] . However, looking into the entropy, we can distin- 
guish them (see Fig.|3li). The universal correction 7 — 1 
is the physical manifestation of the residual coherence 
and can be used as a criterion for judging the coherence 
of the system. 

METHODS SUMMARY 

Here we show Eq. (0 by considering the entropy after 
a given thermodynamic process (see process 1 in Fig. [T]) 
as a function of the timing when the operation starts. 
By invoking the so-called replica trick [|28|l, we have an- 
alytically derived the following relation by taking into 
account the leading-order finite-size correction: 

{S')=S'^^^ + y-l, (3) 

where the average on the left-hand side represents the 
time average on the timing (see Supplementary Infor- 
mation section II). Moreover, we have shown (j/ {S') = 
0{d~^^^) where a is the standard deviation of and 
d is the dimension of the Hilbert space. Since d grows 
exponentially with the number of particles, this means 
that S' = S^g^ + 7 — 1 holds with almost unit probabil- 
ity for any timing of the operation in many-body sys- 
tems. Furthermore, since the finite-size effect 7 — 1 is 
a universal constant, it can be ignored in the thermo- 
dynamic limit independently of the details of the sys- 
tem and processes, and thus we obtain S' = S^^^ > S 
for almost every timing of the operation. It is reason- 
able on physical grounds to consider that the excep- 
tional timings correspond to either before equilibration 
or on revivals after equilibration. Thus, we conclude 
that the entropy never decreases during any thermody- 
namic process between two pure states at equilibrium 
in the thermodynamic limit. 

ACKNOWLEDGEMENT 

Fruitful discussions with Takahiro Sagawa are grate- 
fully acknowledged. This work was supported by 



KAKENHI 22340114, a Grant-in-Aid for Scientific Re- 
search on Innovation Areas "Topological Quantum Phe- 
nomena" (KAKENHI 22103005), a Global COE Program 
"the Physical Sciences Frontier", and the Photon Fron- 
tier Network Program, from MEXT of Japan. T.N.I, ac- 
knowledges the JSPS for financial support (Grant No. 
248408). 



NOTE ADDED 

After completion of this work, we became aware of an 
independent work by Hal Tasaki ll^, where the second 
law of thermodynamics for pure states is addressed in 
terms of the energy. 
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SUPPLEMENTARY INFORMATION 
(Emergent Second Law in Pure Quantum States) 



This supplementary material includes four sections. 
First, we demonstrate how the geometrical concepts uti- 
lized in a qubit system are generalized to higher dimen- 
sions. Second, we show the results presented in Meth- 
ods Summary in the manuscript. Third, we provide 
supplementary numerical data that verify the universal- 
ity of the non-extensive correction to the entropy after a 
thermodynamic process. Last, we apply our method to 
a closely related quantity, the diagonal Renyi entropy 
for which we can also find universal corrections. This 
result suggests that our non-extensive universal correc- 
tion appears for a rather broad class of diagonal en- 
tropies. 



A. GEOMETRICAL INTERPRETATION FOR THE 
HIGHER-DIMENSIONAL HILBERT SPACE 

In this section, we describe the quantum orbit and 
the operation in the high-dimensional "Bloch sphere''. 
Then we illustrate that the entropy never decreases for 
the decohered state. With the results shown in Sec- 
tion |B1 this leads to the entropy increase for pure states. 



f^Pdec 

^dec 




FIG. A.l: Geometry in higher dimensions. In higher dimen- 
sions, the quantum orbit is described in the subspace perpen- 
dicular to the {d — 1) -dimensional plane spanned by the Xa- 
axes {a = 1, 2, . . . , — 1). The operation preserves the distance 

from the origin R = \J Yjf=i ^a/ but decreases the distance 

from the equatorial plane r_\_ = ^Ef=i ^a- The approach to 
the equatorial plane implies an increase in entropy. 



A-1. Bloch sphere 

Let d be the dimension of the Hilbert space. Every 
quantum state is described hy a d x d density matrix p. 
The density matrix p is Hermitian with unit trace. As a 
convenient orthonormal basis of the traceless d x d Her- 
mitian matrices, we adopt the generalized Gell-Mann 
matrices (a = 1, 2, . . . , — 1) [1]. The first d — 1 ma- 
trices are diagonal and analogous to a^, and the others 
are off-diagonal and analogous to Ci and (72/ where CaS 
{a = 1,2,3) are the Pauli matrices. Thus, each density 
matrix p is parametrized by — 1 variables Xa = tr[|OAfl] 
{a = 1, . . . ,d^ — 1) and represented by a point in the 
{d^ — 1) -dimensional Euclidean space. We note that all 
these points, which constitute the "Bloch sphere'' in the 
(d^ — 1) -dimensional Euclidean space, do not form a 
{d^ — 2) -sphere for d > 2 unlike the qubit system for 



A-2. Quantum orbit 

During the time evolution with a fixed Hamilto- 
nian, a quantum state traces a quantum orbit in the 
"Bloch sphere". In the energy eigenbasis representa- 
tion, the first d — 1 components, x^s {a = 1,2, . . .,d — 
1), do not change whereas the others, XaS {a = d,d + 



1, — 1), do. Thus, the quantum orbit is re- 

stricted to a subspace, which is perpendicular to the 
{d — 1) -dimensional plane spanned by the x^-axes {a = 
l,...,d-l) (see Fig.EIJ. 

The decohered state p^ec is the projection of the or- 
bit onto the subspace, in which no coherence exists be- 
tween energy eigenstates. We note that each point on 
the quantum orbit is equidistant from p^^ec because the 
first d — 1 components, XaS (a = 1, 2, . . . , d — 1), and the 

distance from the origin R = ^EfLl"^ are preserved 
under unitary evolutions. 



A-3. Operation 

The operation, which can be implemented by a 
change in control parameters, is represented hy a d x d 
unitary matrix. The change of the energy eigenbasis 
during the operation can also be described by a unitary 
matrix. The combination of these two unitary matrices 
produces a single unitary matrix acting on the quantum 
orbit. Due to the unitarity, the distance from the origin 
R and the distance between pdec arid each state on the 
orbit are preserved. 
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A-4. Monotonicity of the entropy of the decohered state 

The equatorial planeial plane of the Bloch sphere is 
generalized to be the plane where Xa = ior a = 
1,2, . . — 1. The distance from the plane is defined 



by ^ JE: 



=1 -^a- 



We also define r\\ = 



as the distance from the subspace spanned by XaS 



{a = 1,2, . . .,d - 1). Then the relation = r^^ + 
holds. 

The decohered state satisfies ry = by definition and 
thus R = r±. Since R is preserved during the operation, 
r_\_ > r'^ follows. This means that the decohered state 
approaches the equatorial plane during the operation, 
resulting in an increase in entropy. 



B. ENTROPY AFTER A THERMODYNAMIC PROCESS 

In this section, we analyze the entropy after an arbi- 
trary thermodynamic process as a function of the tim- 
ing when the operation starts. The infinite-time aver- 
age involves S^^^ and the universal correction 7 — 1. 
The infinite-time variance is shown to decrease expo- 
nentially with the system size. 



B-1. Setup 



The Hamiltonian is given as follows (see Fig. lB.lj) : 



H{t) 



H 
H' 



op(0 



f < t; 

T < t < T + At; 
T + At < i. 



(B.l) 



where we assume Hop(T) = H and Hop(T + At) = H^ 
Let and {|£^)} be the set of eigenstates of H 

and that of H\ respectively. To ensure equilibration, 
H and H' are assumed to be nonintegrable. The initial 
state is assumed to be pure and thus expanded in terms 
of the energy eigenstates of H: \xp) = En^n l^n)- 
also assume that the number of the superposed eigen- 
states is large and the inverse participation ratio (IPR) 
Q = En I 1^ is of the order oid~^ where d is the dimen- 
sion of the Hilbert space and grows exponentially with 
the system size. The IPR plays a role as the expansion 
parameter in the following discussion. 

At time f = r, each eigenstate has acquired a 
phase proportional to its eigenenergy: |tp(T)) = 
En Cne~^^''^ where Planck's constant is set to unity. 
The time evolution from r to r + At is described by 

a unitary operation V = Texp ^— i J^^ Hop{r + t)dt^, 

where T is the time-ordering operator. Thus the quan- 
tum state at t = r + At is given by \xp{r ^ Ar)) = 
En Cne~^^'''^V \ En). Expressing this state in terms of the 
eigenbasis of H\ we obtain |tp(T + At)) = En^n \^n) 
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FIG. B.l: Control protocol. First, an initial pure state with 
Hamiltonian H undergoes equilibration and reaches equilib- 
rium. Then the operation starts at t = r which smoothly 
transforms H to the final Hamiltonian H\ The entropy af- 
ter the operation, in general, depends on r. However, for r 
greater than the equilibration time, is expected to be inde- 
pendent of T apart from small fluctuations due to a finite-size 
effect. 



where Em ^^'^'^Unm and U, 



{^'m\y\^n) 



Note that UmnS denote the transition amplitudes be- 
tween pairs of the eigenstates before and after the oper- 
ation. 



B-2. Entropy after the operation 

The entropy after the operation depends on the tim- 
ing when the operation is performed, and is given by 



s'(T) = -E 



J \2 



In \c 



f |2 



(B.2) 



Since we are interested in the operation performed on 
the pure state at equilibrium, we calculate S'{t) for suf- 
ficiently large r. To examine this limit, we have only to 
calculate the long-time average on r, or (S^(t))t. Here 

{• • •)t means lim^^oo Jo ' ' ' ^t. We will also show 
below that the fluctuation of S^(t) actually vanishes for 
large systems. To calculate the long-time average of 
each term on the right-hand side of Eq. (|B.2|) , we first 
consider the probability distribution of 



B-3. Moments of the population of each eigenstate after 
the operation 



Here we show that the population of each eigenstate, 
p, after the operation obeys the exponential distribu- 
tion. For this purpose, we analyze the moments of the 



3 



fc-th order: 

{\cT). 



h h 



Y[UnmjU*i.CmjCi, 



(B.3) 



where (pm = £mT (mod 2n) and (• • • ) denotes the av- 
erage over (pm's under the assumptions that each (pm is 
uniformly distributed over [0,27r) and that (pm and (pi 
are independent for m ^ 1. The replacement of (• • • )t 
with (• • • ) is justified for nonintegrable systems, where 
the energy spectrum is sufficiently complex 10]. 
We note the following relation: 



coherence, this state is stationary After the operation, 
it becomes jO^^^ = Vp^Q^^^, and the population of each 
eigenstate |£^) is = jin- Thus, }in corre- 

sponds to the population of for the decohered ini- 
tial state. We note that the factor k\ in Eq. (IB.5D would 
be absent if it were not for coherence between energy 
eigenstates. 

Equation (IB.5I) means that |c^p obeys an exponential 
distribution 



P(\c- 



/ |2^ 



/}ln, 



(B.6) 



E tl^nijl^if) + (corrections). 



(B.4) 



whose mean and variance are }in and respectively 



B-4. Average of the entropy after the operation 

We obtain the long-time average of S' {t) by taking 
the average over the exponential distribution [Eq. (|B.6I) 1: 



where Sy^ is the symmetric group on k elements. It can 
be shown that the correction terms give contributions 
smaller than the first term on the right-hand side of 

Eq. (|B3| roughly by a factor of Qn = Lm {\Unm\^\cm\^)^ • 
Here Qn is a kind of the IPR defined above and 
represents an effective number of the eigenstates that 
can evolve into by the unitary operation V. 

Let us illustrate this for k = 2. Writing out the terms 



POO 

(s^(t))t = -E/ K 

= -^^n In + 7-1 



-die: 



/ |2 



dec 



(B.7) 



explicitly, we have (e-^(^-i+^-2-^^i-^^2)) 



where S^^^ is the entropy of p^g^ and 7 is Euler's con- 
stant. Interestingly, the constant 7 — 1 does not depend 
on details of the system and the operation, whereas S^^^ 



^mim2^mili^m2h 



The third term, whichls the ^^^^ ^dec extensive H, and S < 



/ 

dec 



only correction term in this case, corrects the double 
counting in the first and second term for mi = m2- 
In such a manner, the correction terms in Eq. (|B.4|) 
have more Kronecker's deltas which make the correc- 
tion terms smaller by a factor of Qn- 

We assume that the operation is such that ^ 1. 
This condition is satisfied unless At = 00 or Hop(f) = 
H = H' , because, in large systems, the energy spectrum 
so dense that there should be a number of transitions 
between many-body eigenstates. We note that At = 00 
corresponds to the adiabatic process, which causes no 
transition between many-body eigenstates. Although 
our discussion cannot be applied to this case, we obtain 
S = by definition, which is consistent with the fact 
that the entropy is invariant for quasi-static processes in 
isolated systems. As for Hop(f) = H = H' , the entropy 
is invariant for lack of any operation. 

Substituting Eq. (|B.4|I into Eq. (IB.5II and ignoring the 
correction terms, we obtain 



smce \Umn\^ is doubly-stochastic |H]. Thus, the second 
law is obtained for large systems, where the constant 
7 — 1 can be ignored. 

To make it clear that the universal constant 7 — 1 
originates from the factorial k\ in Eq. (IB.5L or, equiv- 
alently, from the quantum coherence between energy 
eigenstates, we provide another derivation by invoking 
the so-called replica trick |0], which relies on the ana- 
lytic continuation of Eq. (|B.5|) to fc = 0: 



(S(t)).= 



riim 



I{k+2)^l+^-^n 



Y^fln^fin +7-1- 



(B.8) 



(14 



\2k\ 



klpl 



It is clear from this derivation that the universal con- 
stant arises from fc! = r(fc + 1). 



(B.5) 



where }in Lm \ Unm\^\Cm\^- 

To see the physical meaning of finr let us imagine that 
the initial state has no coherence between energy eigen- 
states. The decohered state is represented by a density 
matrix p^ec = En k^p Due to the absence of 



B-5. Variance of the entropy after the operation 

Finally, we justify our conjecture that S'(t) is effec- 
tively independent of t if t is larger than the equilibra- 
tion time and the operation is performed on the pure 



4 



state at equilibrium. For this purpose, we show that system size. 

the variance of S'{t) decreases exponentially with the We begin by calculating the moments for n ^ n, 



\cT\c'nn 



E E 



mi,...,m]^ mi,...,% 17=1 



Y[UnmjU^l.CmjCi. 



1 I ^nfhj^fll^^mj^l^ 



i[LU fx-i+L}^^ 'Pmj-LU 'Pi-q^^ 



(B.9) 



We assume k <k without loss of generality. Substitut- where 2fi(^, ^^,c;z) is the Gauss hypergeometric func- 
ing Eq. (|E1|) into Eq. (|B2| and ignoring the correction tion and Rnn = {vnn?- / [^nUn)- We note that Rnn ^ 

Q <C 1 because En,n knnP = En \cn\^ = Q < 1, while 
YLnfiV-n^n = 1- Thus, the first-order approximation is 
justified, which leads to 



terms, we obtain 

(14 



/ \2k\J \lk\ 



(B.IO) 



where Vnn Em ^nnMfi m\^m? ' are the off-diagonal ele- 
ments of Equation dB.lQI) can be expressed as 



(14 



/ \2k\J \2k\ 



T{k + l)iil T{k + l)iil 2Fi{-k,-k,V,Rnn), 



(14 



/ \2k\J \2k\ 



(B.12) 



(B.ll) gy invoking the replica trick, we have for n n 



(|c'„|2ln|c'„|2|4|2ln|4|2)^ - (|c'„|2ln|c;|2)^(|4|2ln|4|2)^ = J?„, • O ((d"! Ind)^) , (B.13) 



where we assume fin = 0{d ^) since En /^w = 1- Fi" 
nally, we obtain the variance of S\t): 

= E[((KI'in|4l')V-(Kl^in|4l^)? 

+ E[(l4l'ln|4n4l'ln|4l')x 

n,n 

-(|4|^ln|4|^),(|4l'ln|4l')T 

= O • (d In d)^) + O (d^ • d'^ {d In d)^) 

= o(d-^(lnd)^), (B.14) 

where we have used En = 0{d) and Rnn = 0(d~^). 
Since (S'(t))t = O(lnd), the variance vanishes expo- 
nentially with the number of particles N in the system: 

a 



(S'(T))x 



0(d-i/2) = o{e- 



(B.15) 



According to the Chebyshev's inequality, the probabil- 
ity with which S'{t) deviates significantly from (S'(t)) 
becomes exponentially small as the system size in- 
creases ||1.0I- 



C. UNIVERSALITY OF THE NON-EXTENSIVE 
CORRECTION 



We provide further evidence for the universality of 
the correction 7 — 1. The model and the protocol are 
the same as in the manuscript. We focus on the second 
operation and plot S^^ — S^^^^ with the error bar a for 
various values of h (see Fig. IC.ID . Here S^^ and a are 
calculated by the average and the standard deviation 
of S^\r) in the time interval [50,100] in units of h/J. 
We can see that the universal correction coincides with 
7 — 1 independently of the value of h. 
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FIG. C.l: Universality of 7 — 1. S'' averaged over r G 
[50,100]/z// minus S^^^ plotted for various h. Tlie error bar 
is the standard deviation. Note that 7 — 1 is independent of 
h. 



D. THE DIAGONAL RENYI ENTROPY 



v\^here k is an integer larger than 1. The analytic con- 
tinuation k ^ 1 gives the entropy discussed above. We 
also note that S2 = — In Q. By invoking the replica trick 
and the approximation used in calculating the variance 
of the entropy, we obtain 



(S[(T)). = S[ 



k,dec 



InM 
1^' 



(D.2) 



ln(EnF«)/(l 
Renyi-entropy counterpart of 



k) is the diagonal- 

The second term on 
the right-hand side of Eq. (|D.2|I is independent of the 
details of the system and processes and thus universal. 
Taking the limit of fc ^ 1 in Eq. (|D.2|) , we reproduce the 
result (|B7|) . 



Equation (|D.2|) can also be obtained in a different way. 



We introduce the generalized IPR as Q^^^ 
Then, from Eq. (|B.5D , we have 



\2k 



Here we show that the universal correction can also [Q^^^ ' {t))t = k\Q^^^' , (D.3) 

be seen in the diagonal Renyi entropy, which is defined 



as 



In knp^l where Q^^^' dec = En /^n- Taking the logarithm of both 

V ^ ' ^' ) sides and interchange the logarithm and the average on 

1 — k the left-hand side, we obtain Eq. (|D.2|) . 
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